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We use polarised neutron diffraction to study the induced magnetization density of near optimally 
doped Ba(Feo.935Coo.o65)2As2 (Tc=24 K) as a function of magnetic field (1 < fioH < 9 T) and 
temperature (2 < T < 300 K). The T-dependence of the induced moment in the superconducting 
state is consistent with the Yosida function, characteristic of spin-singlet pairing. The induced 
moment is proportional to applied field for fioH < 9 T « fioHc2/6- In addition to the Yosida spin- 
susceptibility, our results reveal a large zero-field contribution M(H — > 0, T — ► 0)/H m 2/3Xnormai 
which does not scale with the field or number of vortices and is most likely due to the van Vleck 
susceptibility. Magnetic structure factors derived from the polarization dependence of 15 Bragg 
reflections were used to make a maximum entropy reconstruction of the induced magnetization 
distribution in real space. The magnetization is confined to the Fe atoms and the measured density 
distribution is in good agreement with LAPW band structure calculations which suggest that the 
relevant bands near the Fermi energy are of the d xz / yz and d xy type. 



I. INTRODUCTION 



gapped structures [4-7, 9-11]. 



The discovery of superconductivity in iron-based mate- 
rials such as LaOi-^F^FeAs [1] has ignited intense inter- 
est in this field. The common features of the iron-based 
superconductors appear to be that they are semi-metals 
which have electron and hole Fermi surface pockets, sep- 
arated by a (tt, 7r) wavevector [2]. Experiments have 
demonstrated the existence of strong spin excitations 
with this same wavevector for superconducting compo- 
sitions. It is widely believed that iron-based supercon- 
ductivity is mediated by these spin excitations resulting 
in superconducting states such as s±, where the sign of 
the gap changes between different Fermi surface sheets. 

Of particular interest in the iron-based superconduc- 
tors is the structure of the superconducting gap and 
the nature of the low-energy electronic states. Pene- 
tration depth [3] and thermal conductivity studies [4] 
of a number of materials (e.g. LaFePO, KFe2As2 and 
BaFe2(Asi_ a: P 3 ;)2) show evidence for low-energy quasi- 
particlc excitations which could be due to nodes in the 
superconducting order parameter. 

Ba(Fei_ a: Co a: )2As2 is a good system to study since 
it is possible to grow large single crystals with ho- 
mogeneous doping. It has been studied widely by 
various probes including angle-resolved photoemission 
spectroscopy (ARPES) [5], scanning tunnel microscopy 
(STM) [6], penetration depth [7, 8], ^SR [9], heat capac- 
ity [10, 11] and thermal conductivity [4]. Even in this 
single system, different gap characters have been pro- 
posed as a function of doping, including fully and nodally 
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Here we use half-polarized neutron diffraction to mea- 
sure the susceptibility and induced magnetization in 
the normal and superconducting states of near optically 
doped Ba(Fei_ x Co x )2As 2 (x = 0.065). Our measure- 
ments shed light on the electronic structure and the na- 
ture of the low energy electronic states in both phases. 
By measuring the flipping ratios of a number of Bragg 
peaks, we are able to extract the spatial Fourier com- 
ponents of the induced magnetization density M(r). In 
a metal, this provides information about the electronic 
states near the Fermi energy. We compare our results 
with a band structure calculation. In addition to mea- 
suring the M(r), we also made a detailed study of the 
temperature and magnetic field dependence of the in- 
duced magnetization by measuring a single Bragg peak 
in detail. 



The paper is organised as follows. In Sec. II we in- 
troduce the polarized beam method used in our exper- 
iment and report our unpolarized structural refinement 
of Ba(Fei_ x Co x )2As2 (x = 0.065). In Sec. Ill we report 
our determination of the real space magnetization den- 
sity M(r) induced by an applied magnetic field. We also 
present the results of a LAPW calculation of the magne- 
tization density distribution. In Sec. IV, we report mea- 
surements of the induced magnetization in the supercon- 
ducting state as a function of magnetic field and tempera- 
ture. We discuss the significance of our observations with 
respect to the superconductivity in Ba(Fei_ a; Co a; )2As2 
(x = 0.065) and other experimental results. This is fol- 
lowed by a summary of our conclusions. 
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II. BACKGROUND 

A. Polarized Neutron Diffraction Studies of the 
Induced Magnetization 

Polarized neutron scattering experiments can directly 
measure the real-space magnetization density M(r) in 
the unit cell, induced by a large magnetic field /iqH . Due 
to the periodic crystal structure, the applied magnetic 
field induces a magnetization density with spatial Fourier 
components M(G), where G are the reciprocal lattice 
vectors, such that 



Atom Position in I4/mmm ~z B (A 2 ) 
Ba 2a (0 0) 0.023(16) 
Fe/Co 4d (| i) 0.110(10) 
As 4e (0 z) 0.35352 083(11) 

a = 3.952(2) A,c = 12.911(15) A 
g=22.4(1.2) rad' 1 



FABLE I. Structural parameters of Ba(Fei_ a; Co a: )2As2 (x = 
0.065). Parameters are obtained from least-squares refine- 
ment of integrated intensities measured at T — 30 K on D9. 
The B-factor is related to the mean squared atomic displace- 
ment (u 2 ) by B = 8ir 2 (u 2 ). g is the width parameter of the 
mosaic distribution [16, 17]. 



M(r) = —^2 M(G)exp(-iG • r) 
Vq g 



(1) 



and z/q is the volume of the unit cell. The Fourier com- 
ponents of the magnetization density are given by 

M(G) = f M(r) cxp(iG • r)dr. (2) 

Amit cell 

Neutrons interact with matter through the strong nu- 
clear force and electromagnetic interaction. For neutrons 
with initial and final spin polarisations <7j and cry, the to- 
tal scattering cross section is 



(da 
\dtt 



| ir 

2fi B 



G x |M(G) x G| 



F N (G)\a f )f 



(3) 



where -fr — 5.36 x 10 -15 m and Fjy(G) is the nuclear 
structure factor. The sign of the first (magnetic) term in 
Eq. 3 can be changed by reversing the incident neutron 
polarization. Thus we are able to isolate the interference 
term between the nuclear and magnetic scattering. In 
this experiment we measure the flipping ratio R, defined 
as the ratio of the cross-sections with neutrons parallel 
and anti-parallel to the applied magnetic field. Because 
the induced moment is small, the experiment is carried 
out in the limit (yr /2fj, B )M(G)/F N (G) < 1. In this 
limit, the flipping ratio derived from Eq. 3 is, 



R 



\F N (G) - ( 1 r /2 f i B )M(G)\ 2 _ i 2 7 r M(G) 
|^(G) + ( 7 ro/2/i B )M(G)| 2 ~ 



Ms F N (G)' 

(4) 

As the nuclear structure factors Fjv(G) are known from 
the crystal structure, Eq. 4 directly gives M(G). 



B. Experimental Details 

The Ba(Fei_ :c Co a ;)2As2 (x = 0.065) single crystal used 
in this study was prepared by a self- flux method [12], had 
approximate dimensions 6 x 1.5 x 0.2 mm 3 and a mass of 
~1.8 mg. Similar samples were used in our previous stud- 
ies [13, 14]. Resistivity and magnetization measurements 



on crystals from the same batch identified the super- 
conducting transition temperature T c (onset)=24 K and 
showed no evidence of magnetic order down to 2 K. The 
bulk susceptibility in a 5 T field measured using a SQUID 
magnetometer was \ab = 1-22 x 10~ 3 Mb T _1 f.u. -1 . We 
note that the expected upper critical field applied in the 
ab plane for this composition is H C 2, a b =55 T [15]. 

Neutron scattering experiments were performed at the 
Institut Laue-Langevin, Grenoble, France. An initial un- 
polarized structural refinement was performed using the 
4-circle D9 spectrometer with A = 0.837(1) A. The re- 
sults are shown in Table I. Polarized beam measurements 
of the flipping ratio were made on the D3 spectrometer. 
The sample was mounted on a thin aluminium post with 
the [110] direction vertical and parallel to the applied 
field. Data was collected with an incident wavelength 
A = 0.825 A and a 0.5 mm Er filter (to reduce higher or- 
der contamination in the incident beam) . Flipping ratios 
for equivalent reflections were collected, averaged and 
corrected for the finite beam polarization and extinction 
effects. The sample was cooled through T c at each field 
measured when collecting data in the superconducting 
state. 



III. INDUCED MAGNETIZATION 
DISTRIBUTION 

A. Results 



We measured the magnetization induced by a 9 T mag- 
netic field applied along the [110] direction in the normal 
state at T=30 K. Table II shows the measured flipping 
ratios under these conditions and the extracted Fourier 
components of the magnetization density. Fig. 1 shows 
|M(G)| plotted against sintf/A = |G|/4tt. If Fe were 
the only magnetic atom, then the sign of M(G) is de- 
termined by the geometric structure factor for the Fe 
atoms and |M(G)| is the effective form factor of the Fe 
atom. The solid line in Fig. 1 is the standard calculated 
isotropic atomic form factor for Fe 2+ [19]. Deviations 
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TABLE II. For each Bragg reflection (hkl), the table shows: 
sin 8/\, the measured flipping ratio R, structure factor Fn (G) 
calculated from the structure in Table I, and the determined 
magnetic structure factor M(G) [18]. Data collected at T = 
30 K and fi H = 9 T. 



Ba(Fe 1x Co x )As 2 (x=0.065) \l q H=9 T 
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FIG. 2. (left) The conventional tetragonal unit cell of 
Ba(Fe,Co)2 As2. We measure the magnetization density in- 
tegrated perpendicular to the (llO)-type plane shown. The 
magnetization density shown is the result of the VCA cal- 
culation shown in Fig. 4. (right) Projected magnetization 
distribution reconstructed from the experimental data in Ta- 
ble. II. The reconstructed magnetization map shows j of the 
area of the plane shown in the left panel. 



FIG. 1. The magnetic structure factors |M(G)| measured for 
a 9 T field applied along the [110] direction. The solid line is 
the Fe 2+ form factor [19] scaled to the measured susceptibility. 



from an isotropic form factor are expected at larger 9 or 
G|. Our results are in broad agreement with a recent 
study of Ba(Fei_ K Co K )2As 2 (x = 0.066) [20]. The main 
differences between the present data and the data pre- 
sented in Ref. 20 are: (i) the present study has higher 
statistical accuracy (ii) the value of M(G) for G=(002) 
is larger relative to the other G points in the present 
study. Our extinction model shows that the (002) peak 
has the largest extinction correction factor of 1.26. 

We used the maximum entropy method (MaxEnt) [21- 
24] to make a model-free reconstruction of the magnetiza- 



tion density in real space. Flipping ratios for reflections 
of the type (hhl) were collected, this allowed the recon- 
struction of the magnetization density projected down 
the [110] direction onto the (110) plane as illustrated in 
the left panel of Fig. 2. The result of the reconstruc- 
tion is shown in the right panel of Fig. 2. As expected, 
the magnetization density is localized mostly on the Fe 
atoms. The magnetization "cloud" appears to be slightly 
extended along the (110) directions. Our results are in 
broad agreement with Ref. [20], however we observe no 
significant magnetization density on the Ba sites. A re- 
cent study of the paramagnetic state of the parent an- 
tiferromagnet BaFc2As2 [24] is also broadly consistent 
with our results. The main difference in the BaFe2As2 
case is that the magnetization extends more towards the 
As atoms. 
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B. Electronic Structure Calculations 



Induced form factor measurements have been widely 
used to determine the nature of the electrons responsi- 
ble for paramagnetism in solids. In metals, the induced 
magnetization arises from a redistribution of electrons 
between up and down states near the Fermi energy. Thus, 
we probe the nature of the electronic wavefunctions for 
states near the Fermi energy. In order to understand our 
results further, we have carried out electronic structure 
calculations using the WEIN2k package [25]. 

The linear augmented plane-wave (LAPW) method 
[26] was used to obtain the electronic structure and spin 
density of Ba(Fei_ :E Co :E )2As2 (x = 0.065). We used a 
full-potential LAPW method with the generalized gradi- 
ent approximation (GGA). In the case of doped compo- 
sitions, we used the virtual crystal approximation (VCA) 
[27]. We used the lattice parameters and atom posi- 
tion for the Ba(Fei_ x Co a: )2As2 (x = 0.065) structure 
shown in Table. II for both compositions. The muffin-tin 
radii were chosen to be 2.5, 2.37 and 2.11 atomic units 
(a.u.) for Ba, Fe and As respectively, with the quantity 
i?MT^max set to 7, where Rmt is the smallest muffin-tin 
radius and if max is the plane wave cutoff. For integra- 
tions we used 726 fc-points in the irreducible Brillouin 
zone. In the LAPW method, the charge (spin) density 
is represented by a plane wave expansion in the inter- 
stitial region and as a combination of radial functions 
times spherical harmonics inside the mufhn-tin spheres. 
Fig. 3 shows the bandstructure calculated in the tetrag- 
onal phase for Ba(Fei_ a; Co z: )2As2 for x=0 and a;=0.065. 
Our calculations generally agree with others in the lit- 
erature [28-31], in particular, they show that the states 
near the Fermi energy are predominately of d XZtVZ and 
d xy character. 

Spin polarized calculations were carried out in the 
tetragonal state with the ferromagnetic magnetization 
constrained to be 0.01 \ib per unit cell to mimic the ef- 
fect of an applied magnetic field. The results were then 
scaled to the measured magnetization at fioH =9 T for 
comparison purposes. The results are shown in Fig. 4. 
The x = 0.065 calculation shows reasonable agreement 
with our maximum entropy reconstruction in that the 
magnetization density is extended along the (110) direc- 
tion. However, the reconstruction from our experimen- 
tal data does not show the two maxima along (£,£, 1/2) 
predicted in the calculation. This is presumably be- 
cause our data have insufficient Fourier components to 
resolve these features. We also carried out a calculation 
for BaFe2As2 (x = 0) in the paramagnetic tetragonal 
state with the same structural parameters (but no Co 
potential) to demonstrate the sensitivity of the magne- 
tization distribution to the bandstructure. Notice that 
the change in electronic structure between x = and 
x = 0.065 causes a rotation of the calculated pattern in 
Fig. 4. 



IV. TEMPERATURE AND FIELD 
DEPENDENCE OF THE MAGNETIZATION 

The measurement of the field and temperature de- 
pendence of bulk magnetization M(H,T) in the mixed 
state of a superconductor provides information about the 
nature of the superconductive pairing. Thermal con- 
ductivity k(H,T,9) [4] and specific heat measurements 
C = "j{H)T [32] in the mixed state have been extremely 
useful in probing the low energy quasi particles and iden- 
tified the gap structure of a number of superconductors. 
In particular, the field dependence of the electronic con- 
tributions to k and C in the T — >• limit is sensitive to 
the symmetry of the superconducting gap function [33] 
(see Fig. 5). Complementary information is contained 
in M(H,T). However, studies of the bulk magnetization 
in the mixed state of superconductors are not possible 
by conventional means, e.g. SQUID magnctometery, be- 
cause of the presence of a large diamagnctic contribution. 
Polarized neutron diffraction and NMR Knight shift mea- 
surements [35, 36] can be used to make magnetization 
measurements in the mixed state. The polarized neu- 
tron diffraction technique used here is unique because it 
directly measures the total magnetization including spin 
and orbital contributions. This technique was first used 
by Shull and Wedgwood in 1963 to study V^Si [37] and 
has subsequently been applied to such varied supercon- 
ductors as UPt 3 [38], YBa 2 Cu30 7 _ (5 [39] and Sr 2 Ru0 4 
[40]. 



A. Results 

We used the (002) Bragg peak to make our measure- 
ment since this requires little correction to give the G = 
magnetization, M(G = (002))/M(G = 0) = 0.979 for an 
isotropic Fe 2+ ion. Fig. 6 shows the temperature depen- 
dence the induced magnetization of Ba(Fci_ 2; Co 2 ;)2As2 
(x=0.065) for hqH = 9 T. The temperature dependence 
shows the characteristic Yosida dependence expected for 
singlet pairing [41], this is consistent with measurements 
of the NMR Knight shift in the same compound [35, 36]. 
There is a large susceptibility in the T — > limit. A large 
residual contribution is also observed in V 3 Si [37] and this 
has been attributed to the van Vleck contribution [42]. 
The presence of large van Vleck contribution has also 
been inferred [43, 44] from NMR knight shift measure- 
ments on other iron-based superconductors. Fig. 7(top 
panel) shows the field dependence of the induced magne- 
tization. M(G) oc H over the field range of the present 
experiment ((j,qH < 9 T). When we plot the susceptibil- 
ity M(G)/H [see Fig. 7(bottom panel)], we find that the 
value as T — > and H — > is about | of the normal 
state value measured at T c in the present experiment. 
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FIG. 3. (Color online) (bottom panels) The linear augmented plane-wave (LAPW) bandstructure calculations for 
Ba(Fei_ a; Co :E )2As2 for x=0 (left) and a;=0.065 (right) in the virtual crystal approximation (VCA). Colors indicate the dominant 
orbital of the k state, (top panels) Calculated electronic density of states. 



B. Interpretation 

We measured the static response at finite wavevector 
G to an (approximately) uniform magnetic field, that 
is, X (G) = X (G,0) = M(G)/ff(G = 0) [18, 45]. The 
atomic nature of solids means that the induced magneti- 
zation M(r) is spatially modulated on an atomic scale 
(see Fig. 2). Neutrons diffract from this modulation. 
In the mixed state of a superconductor, there is an ad- 
ditional diamagnetic magnetization (which gives rise to 
the vortex lattice), which is not detected in the present 
experiment. The signal from the vortex lattice is only 
present at small wavevectors (scattering angles) and can 
be studied by neutrons using small angle scattering tech- 
niques. In the following discussion wc do not include the 
superconducting diamagnetic response. 

The magnetic susceptibility in <i-band metals has sev- 
eral components: atomic diamagnetic, van Vleck ( "or- 
bital" or "interband") and spin. The atomic diamagnetic 



contribution in BaFe2As2 has been estimated [24] to be 
small and is neglected here. Only the spin contribution 
is expected to be suppressed by singlet Cooper pairing, 
thus we write the spatially averaged induced magnetiza- 
tion in the mixed state as: 

M = X oTbH + M spin (T,H), (5) 

where %orb is the orbital susceptibility and M sp i n is the 
spin magnetization. In an s-wave superconductor the 
density of states due to the introduction of vortices is 
oc NpH/H C 2 [46, 47]. Thus the spin magnetization 
should vary as M spin (T -» 0,H) oc H 2 . The tem- 
perature dependence is given by the Yosida function 
[41], M S pi n (T, H -> 0) cx Y(T). The present experi- 
ments (Fig. 7) show that the differential susceptibility 
dM/dH has a large finite value in the H — > limit and 
dM/dH ps constant for /j,qH < 9 T. This is consistent 
with the finite x(T — > 0,H — > 0) response being due 
to a van Vleck contribution. It should be noted that 
there is also NMR evidence for a residual susceptibility 
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FIG. 4. Magnetization density of Ba(Fei_ ;E Co :E )2As2 calcu- 
lated using the LAPW method and VCA approximation for 
x — (left panel) and x — 0.065 (right panel). Calculations 
were carried out with a small fixed ferromagnetic moment to 
mimic the effect of an applied field. 




FIG. 5. Schematic field dependence of the low temperature 
magnetization M, specific heat C, and thermal conductivity 
k for an s-wave superconductor. The behavior for C and k 
are based on VsSi [34]. 
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FIG. 6. The temperature dependence of the suscepti- 
bility and induced moment of VsSi [37] (top panel) and 
Ba(Fei_ x Co :E )2As2 1=0.065 (bottom panel) measured using 
polarised neutron scattering. The solid lines are the Yosida 
behaviour expected for a singlet order parameter. 



in BaFe 2 (As .67Po.33)2 [48] and BaFc(Fc .9 3 Coo.o7)2As2 
[36]. Specific heat measurements also suggest that there 
can be a sizeable residual quasiparticle density of states 
in BaFe(Fei_ l: Co l: )2As2 [10, 11, 49]. Wc cannot rule out 
the possibility that this is related to the residual suscep- 
tibility observed by neutron scattering. 

There has been considerable debate about the nature 
of the superconducting gap in Ba(Fei_ a; Co a ;)2As2. In 
principle, a detailed measurement of the T-dcpendence of 
the induced moment could be used to distinguish between 
different models for the gap. Unfortunately, the statisti- 
cal noise in the present data is relatively high. Thus we 
make only a basic comparison with a singlet s-wave state. 
Within the statistical error of our data, the temperature 
dependent component of the induced moment Af sp j n is 
well described by a Yosida temperature dependence (see 
Fig. 6) with A = 1.78/tflT = 41 K= 3.5 mcV. 
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V. CONCLUSIONS 



Wc have used a polarised neutron diffraction technique 
to measure the induced magnetization density of near 
optimally doped Ba(Feo.935Coo.o6s)2As2 (1^=24 K) as 
a function of magnetic field and temperature. The in- 
duced magnetization is confined to the Fe atoms with an 
oblate distribution spread out in the a — b plane. The 
distribution is in reasonable agreement with a full po- 
tential LAPW band structure calculation which suggests 
that the relevant bands near the Fermi energy are of the 
dxz/yz and d xy type. 

The T-dcpcndcncc of the induced moment in the su- 
perconducting state is consistent with the Yosida func- 
tion characteristic of spin-singlet pairing, and the in- 
duced moment is proportional to applied field for [IqH < 
9 T k /i -ff C 2/6. We observe a large residual suscepti- 
bility M(H -> 0,T 0)/H k 2/3xnormai- This is most 
easily interpreted as being due to the van Vleck contribu- 
tion present in other <i-band systems, but may also signal 
a residual quasiparticle density of states. 
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